MeV Dark Matter: Has It Been Detected? 
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We discuss the possibility that the recent detection of 511 keV 7-rays from the galactic bulge, as 
observed by INTEGRAL, is a consequence of low mass (~MeV) particle dark matter annihilations. 
We discuss the type of halo profile favored by the observations as well as the size of the annihilation 
cross section needed to account for the signal. We find that such a scenario is consistent with the 
observed dark matter relic density and other constraints from astrophysics and particle physics. 
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I. INTRODUCTION 

Recently, observations of a bright 511 keV 7-ray line 
from the galactic bulge have been reported [1]. This de- 
tection, made by the SPI spectrometer on the INTE- 
GRAL (INTErnational Gamma-Ray Astrophysics Labo- 
ratory) satellite, suggests a spherically symmetric distri- 
bution of 511 keV 7-rays with a full width half maximum 
of ~ 9° (6° - 18° at 2a confidence). The flux of the 511 
keV line emission from this region has been measured to 
be 9.9^21 X 10"'* phcm^^ s^^ with a width of about 3 
keV, in good agreement with previous measurements [2] . 

The source of positrons in the galaxy is the subject 
of a great deal of debate. Some suggested sources in- 
clude neutron stars or black holes [3], radioactive nuclei 
from supernovae, novae, red giants or Wolf-Rayet stars 
[4] , cosmic ray interactions with the interstellar medium 
[5], pulsars [6] and stellar flares. Despite this considerable 
list of potential sources, the origin of these positrons is 
still unknown today. Indeed, the most generous sources 
seem to be type la supernovae, however the fraction of 
positrons that escape from such an event, depending on 
internal mixing immediately after explosion and on the 
magnetic field strength and topology, is subject to debate 
[7]. 

Massive Wolf-Rayet stars (hypernova), of the SN2003 
type [8], exploding in the galactic center could do the 
job [9], but their rate is unknown. Also, it is not known 
for certain that the positrons injected could fill the whole 
galactic bulge, even if a bipolar galactic wind is produced 
by star bursts [10]. Indeed, even if a "galactic positron 
fountain" develops [11], the annihilation rate at high al- 
titude is too low, due to the small density of the wind, 
to explain the extension of the 511 keV source [12]. 

Observations of the cosmic microwave background, the 
primordial abundances of light elements and large scale 
structure have revealed that a great deal of the mass of 
our universe consists of dark matter (DM) [13]. Despite 
this large and growing body of evidence, we are still ig- 
norant of the nature of DM. 

The most commonly discussed DM candidates are 



weakly interacting particles with masses typically in the 
range of a few GeV to several TeV, the most popular 
example of this being the lightest supersymmetric parti- 
cle [14] made stable by the virtue of R-parity [15]. This 
particle, along with other fermionic DM particles, is not 
expected to be significantly lighter than a few GeV if 
they are to provide the measured density of DM in the 
universe [16]. This constraint does not strictly apply to 
all particle candidates, however. Light particles (~MeV- 
GeV), in some cases, may be able to satisfy relic density 
and other constraints [17,18]. In particular, light scalar 
DM candidates are an interesting and viable possibility 
[18]. 

In this Letter, we argue that light DM particles ('--^l- 
100 MeV) annihilating into e"'"e~ pairs in the galactic 
bulge may be the source of the observed 511 keV emission 
line. 



II. POSITRON PROPAGATION AND 
ANNIHILATION 

If DM particles of ~1-100 MeV mass annihilate in the 
galactic bulge into electron-positron pairs, the resulting 
positrons then travel with their propagation being dom- 
inated by ionisation losses. This energy loss rate is ap- 
proximately given by [19] 
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where T is the positron's Lorentz factor and Nh is the 
number density of target atoms in the galactic bulge. 
This rate yields a stopping distance of ~ 10^'' and ~ 10^^ 
cm for positrons of MeV and 100 MeV energy, respec- 
tively. For microgauss scale magnetic fields, the Larmor 
radius is on the order of 10^^ or 10^ cm for energies of 
100 MeV and 1 MeV, respectively. Considering a sim- 
ple random walk, the positrons are roughly confined to a 
distance of ~ \J Rsto-p x -RLarm, which is about a parsec, 
or less. Therefore the positrons are easily confined to the 
galactic bulge as they lose their energy. 
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The cross section for the pair annihilation of a positron 
with an electron at rest is simply, CTpairann — (lbn)/r. 
This leads to a positron mean free path of 
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10^5 cm, 
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where is the electron density. Thus, we see that a 
positron's stopping distance is typically much shorter 
than its mean free path. Annihilations are, therefore, 
expected to occur primarily for thermalized positrons, 
producing a 511 keV line. For particles of ^100 McV 
mass, it may be possible that a substantial fraction 
of the positrons could annihilate before reaching non- 
relativistic energies. For this reason, lighter DM particles 
(~1-10 MeV) may be somewhat favored. If the electron 
temperature in the galactic bulge is too low, positronium 
formation may dominate, and results in a narrow line 
(25% of the time) or 3-photon continuum (75% of the 
time), depending on the spin state of the positronium. 
The data sugcsts that annihilations through positronium 
dominate, resulting in a narrow 511 keV line [20]. 



III. THE DARK MATTER HALO PROFILE 

Assuming an equal number of DM particles and anti- 
particles, the DM annihilation rate is proportional to the 
DM density squared. Independent of the particle mass 
or annihilation cross section, the distribution of 511 kcV 
7-rays from the galactic bulge can be used to constrain 
the DM halo profile of our galaxy. 

The halo profile is typically parameterized by 
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where a, (3 and 7 are unitless parameters set by the choice 
of halo model and a is the distance from the galactic 
center at which the power law breaks, also set by the 
halo model. For the inner kiloparsec or so, the region we 
are concerned with, r is much smaller than a, and the 
profile description simplifies to 



p(r) cx — 



(4) 



The flux of annihilation products observed in a region of 
the sky is often described by the function J(^') which 
is proportional to the local DM density squared, inte- 
grated over the line of sight in the direction For such 
a function to be meaningful, it should be averaged over 
the angular resolution of the experiment. 

The distribution of events seen by INTEGRAL ap- 
pears to fit reasonably well to a gaussian with full width 
half maximum of ~ 9°, with a 2a confidence interval 
of 6° — 18°. In figure 1, we compared distributions of 
J(^') averaged over the 2° angular resolution of the SPI 



spectrometer, for a variety of halo profiles. Keeping in 
mind that the angular resolution used here is approxi- 
mate and that there are large error bars for the width 
of the observed 7-ray distribution, this figure indicates 
that halo profiles with the modest slope of 7 ~ 0.4 — 0.8 
provide the best fit to the observed distribution. Such 
a value corresponds to a mild cusp less steep than the 
NFW distribution based on N-body simulations [21], but 
more steep than the model proposed by Kravtsov et al. 
which is based on observations of dwarf and low surface 
brightness galaxies [22]. 
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FIG. 1. The angular distribution of 7-rays from DM an- 
nihilation averaged over the 2° angular resolution of the SPI 
spectrometer on INTEGRAL for several halo profiles. SPI's 
observation indicates a full width half maximum of 9° with a 
6° — 18° 2(7 confidence interval. Shown as vertical dashed and 
dotted lines are the central value and 2a limits of the angular 
widths found by SPI. To agree with this data, a halo model 
with 7 ~ 0.4-0.8 is favored. 



IV. THE ANNIHILATION CROSS SECTION 

If DM is made up of particles with masses below the 
muon and pion masses, low velocity annihilations can 
only produce electron-positron pairs, or perhaps photons 
or neutrinos. Assuming the former dominates, for each 
pair of DM particles which annihilate, one positron is 
produced, which eventually annihilates with an electron 
producing a pair of 7-rays. The flux of 7-rays produced 
in this scenario is given by 
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pb 
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where mdm is the mass of the DM particle, av is the 
annihilation cross section multiplied by the relative ve- 
locity in units of c, Afi is the solid angle observed and 
J(ArJ) is the average of J{'^) over that solid angle. We 
consider a solid angle of 0.02 sr, corresponding to a 9° 
diameter circle. Again, we average over the 2° angular 



2 



resolution of the SPI spectrometer. For this sohd angle 
and resolution, for an NFW halo profile (7=1), we get 
J(0.02sr) ~ 187. As we decrease the profile's slope, we 
get values of 7(0.02 sr) equal to 82.1, 37.3 and 17.1 for 
7=0.8, 0.6 and 0.4, respectively. 

Setting the calculated flux equal to the observed flux 
of 9.9 X 10~^phcin~^ s~^, we arrive at 

f^Vi^Vj(0.02sr)- 0.003. (6) 

Using the width of the observed distri- 
bution (J(0.02sr) ~ 10-100), we can now estimate the 
required annihilation cross section. For a 1 MeV mass, 
av must be ~ 10"'' to 10~^ pb for annihilations in the 
galactic bulge. For 100 MeV mass, roughly 0.1 to 1 pb is 
required. 

V. RELIC DENSITY AND ASTROPHYSICAL 
CONSTRAINTS 

A. Relic Density 

We will now attempt to address the question of: is the 
range of annihilation cross sections found in section IV 
consistent with constraints from relic density and other 
measurements? 

Assuming no coannihilations, the annihilation cross 
section required for a thermal relic is approximately 

where XF=mi^/Tp ~ 17.2-|-ln(5f/5j,)-|-ln(mDM/GeV)-|- 
Iny^xp ^ 12 — 19 for particles in the McV-GcV range, 
and g and 51* are the number of internal and relativistic 
degrees of freedom, respectively. 

B. Astrophysical Constraints 

DM annihilations which occur during the recombina- 
tion epoch can yield substantial fluxes of highly red- 
shifted 7-rays, which may be possible to be observed to- 
day. It has been shown [17], however, that for the class 
of particle we consider in this letter, such constraints are 
not violated. 

Another possible constraint comes from our knowledge 
of nucleosynthesis. A DM candidate may yield too much 
'^He photodissociation if very large S-wavc cross sections 
are present [23]. However, we have checked this con- 
straint and found that the cross sections we consider are 
consistent with all light element abundances. 



C. Assessment 

Combining the constraints from the INTEGRAL mea- 
surement and the relic density calculation, we can assess 
the prospect of a light DM particle being responsible for 
the 511 keV emission in the galactic bulge. The anni- 
hilation cross section needed for -^McV DM to provide 
the correct relic density is somewhat larger than what 
the INTEGRAL data appears to require. This apparent 
discrepancy can be solved, however, if aw oc w^. We can 
write the cross section as av ~ a + bv"^, where a and b 
are some constants related to the S and P-wave terms. 
If a 3> &, then the annihilation cross section is roughly 
velocity independent and the constraints discussed above 
can only be reconciled in the heavy portion of the range 
we consider 100 MeV). If the cross section is S-wavc 
suppressed {b :» a), however, considerably lighter DM 
particles 1 MeV) are possible. 

VI. POSSIBLE MODELS AND PARTICLE 
PHYSICS CONSTRAINTS 

Although DM is typically assumed to be heavier, 
masses below 1 GeV can satisfy eq. 7, despite the Lce- 
Weinberg limit, provided that their annihilation cross 
section is almost independent of the DM mass or if an- 
nihilations involve another light particle with weak cou- 
plings to ordinary matter. The DM candidate itself must 
also have suppressed couplings to the Z. 

The first possibility is easy to achieve with scalar 
DM particles annihilating through the exchange of heavy 
particles (>100 GeV, to evade LEP constraints), non- 
chirally coupled to ordinary fermions. In this case the 
a-term in the annihilation cross section is similar in size 
to the 6-term and the annihilation cross section is only 
weakly dependent on velocity. This scenario can provide 
the observed rate of 7-rays and satisfy eq. 7 for masses 
^100 MeV, but not for considerably lighter particles, 
which require S-wave suppressed cross sections. 

The second possibility relies on the exchange of neutral 
and light particles. Assuming, for simplicity, that the 
DM is a scalar, this exchanged particle could be a new 
fermion with negligible or even null couplings to the Z. 
Alternatively, it could also be a new gauge boson from an 
extra C/(l), for example, called a U boson. Light, neutral 
fermions provide a ^-independent cross section but the 
final state of the annihilations (at tree level) would be 
neutrinos rather than e+e" or 7-rays. U boson exchange, 
on the other hand, provides an annihilation cross section 
into e+e~ pairs, which is S-wave suppressed (b » a). 

The U boson properties can be constrained through its 
additional contributions to the muon or electron anoma- 
lous magnetic moments. One finds that light DM par- 
ticles can provide the correct relic density by selecting 
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appropriate values of the coupling tZ-dm-dm* once the 
coupling U-f-f is determined from g — 2 [18]. 

From eq. 7, av ^ pb at the freeze-out epoch for 
mdm ^MeV (xp ^ 12 and ^ 10). Given the de- 
composition: av = a + bv^ (a = for a C/), we find 
6 ~ 10 pb, thus bv'^ ~ pb at the freeze-out epoch (since 
v"^ = Tfo/nidm ^ 1/12), and therefore, av ^ 10~^pb in 
our galaxy. For U boson exchange, therefore, only very 
light DM particles (-^MeV) are viable DM candidates. 

Both the U boson and the DM particle would evade 
colliders limits despite their small masses, as studied in 
[18]. Indeed, the direct production of U bosons (e+e~ 
•yU) is likely to remain imobserved due to the large 
background associated with the photon pair production. 
Also, the cross section e+e~ 7 dm dm* appears to 
be~6xl0-« i^)' (J^)-^pb while the 

e+e~ ^ yvV cross section is ~ 10^'^ ( cfv) P^' H^i^ce, 
the existence of light DM could have escaped low en- 
ergy experiments searching for anomalous single photon 
events (especially if the cross section for e+e^ dm dm* 
is much smaller than the difference between the theoret- 
ical prediction of the cross section e^e~ vv for 3 neu- 
trino families, and its measurement at the Z pole) [24]. 

The two-body decay of the U boson could provide a de- 
tectable signal, but given our couplings, the main channel 
is likely to be into two DM particles (in which case the 
decay would be dominated by invisible modes). On the 
other hand, a [/-boson could, perhaps, be observed in 
nuclear transitions [25]. 

We also note that, in addition to positronium forma- 
tion, higher order annihilation diagrams may contribute 
to a continuum component of the 7-ray spectrum. 

VII. CONCLUSIONS 

We find that the recent observations of the SPI spec- 
trometer on INTEGRAL of 511 keV 7-rays from the 
galactic bulge are consistent with the annihilation of light 
(~1-100 MeV) dark matter particles. The angular width 
of the observed emission favors a mild cusp in the galactic 
halo profile with 7 ^ 0.4-0.8. 

For ~MeV dark matter particles, an S-wave sup- 
pressed annihilation cross section is required. Such a 
cross section would be predicted for interactions in which 
a new, light gauge boson were exchanged. For heavier 
dark matter particles (~100 MeV), a nearly constant 
cross section is needed. For this case, the exchange of 
new, heavy fermions could be responsible. 
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